Field-induced p-wave Superconductivity in Mesoscopic Systems 
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By using Bogoliubov-de Gennes equations, we study superconducting (SC) states in a quasi 2-dimensional 
system of radius R. It is shown that no vortices exist in s-wave SC samples with R < i? c ~ £(0), the T = 
coherence length. We predict that chiral p-wave states exhibit superconductivity for R < R c only in the presence 
of a vortex with opposite chirality. This induced SC phase is a consequence of non-zero chirality of the pairing 
order parameter and implies the presence of chiral edge currents. Our study may be applied to sharply probing 
the pairing symmetry of unconventional superconductors. 
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One of the most fascinating areas in condensed matter 
physics is the study of unconventional superconductors that 
spontaneously break time reversal symmetry [1], such as p+ip 
and d + id superconductors. Due to their non-trivial topology, 
such systems may support exotic objects such as half-quantum 
vortices and zero energy modes of Majorana quasi-particles 
[2]. A key challenge is to identify the pairing symmetry of 
such superconductors. Currently, the Kerr effect and muon 
spin relaxation experiments have played an important role in 
detecting broken time-reversal symmetry. 

Here we propose that mesoscopic systems may provide 
novel platforms for the study of unconventional superconduc- 
tivity due to the experimental availability and intriguing prop- 
erties in small length scale. Recent advancements in nano- 
fabrication have made it possible to study superconductivity 
in mesoscopic samples with size comparable to the coher- 
ence length £ [3]. As the dimension of a superconductor is 
reduced to order of the SC coherence length, the effect due to 
the edge and the vortex core becomes important. The former 
enforces vanishing quasi-particle amplitudes, while the latter 
requires a topological phase factor on the pairing potential, 
whose amplitude vanishes at the center of the vortex core. All 
these have provided motivation for the research of few-vortex 
physics [4]. 

In this letter, we use BdG equations to study SC states in a 
disk of radius R, comparable to the coherence length £. This 
microscopic theory allows to study the system at zero tem- 
perature and with size smaller than £ quantitatively, and also 
enables us to evaluate observables in detail. Our results are 
further supported by intuitive analysis based on the (Ginzberg- 
Landau) GL theory. To be consistent with previous studies 
based on the GL theory [5], we begin with s-wave supercon- 
ductors, and confirm that the vortex state can only be gener- 
ated with system size larger than a critical value R c . For the 
unconventional paring, we predict that chiral p-wave states ex- 
hibit superconductivity for R < R c only in the presence of 
a vortex with opposite chirality, namely an induced SC phe- 
nomenon. A half-quantum vortex in the equal spin pairing 
p-wave state is also studied. The relevance to the possible p- 



wave superconductor Sr2Ru04 is discussed. 

We start with the s-wave case. The BdG equations read 
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In Eq.l, Ui(r) and u,(r) form the two-component wavefunc- 
tion of quasi-particles corresponding to energy Ei in the 
SC state. A(r) is the pairing potential satisfying the self- 
consistent equation 



A(r)= 9 J2<r)v*(r)[l-2f(E i )}, 

Ei<\ 
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with g the coupling constant, f(Ei) the Fermi distribution 
function, and A an energy cut-off. ho(r) is the single elec- 
tron Hamiltonian 
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where m is the electron mass, /i is the chemical potential. 
A(r) is the vector potential determined by the Maxwell's 
equation V x V x A = —j, with the supercurrent density 
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We work in a polar coordinate system (r, &), and consider a 
natural boundary condition Ui(r = R,6) = u,(r =R,6) = 0. 
Note that we neglect the z-dependence on u, v, and A due 
to the fact that effective mass along z-direction is large and 
the thickness of the system is small enough [6]. The single 
particle wavefunction of hg (with A = 0) corresponding to 
the eigenvalue eji is 
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where I is the angular momentum, J/ (x) is the I th order Bessel 
function of the first kind, and otji is the j th zero of J/ (x) [7]. In 
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the presence of pairing, we have Ui(r) = Ui(r)e lW and Vi(r) = 
Vi(r)e l ( [ ~ n ^ 6 , with n the vorticity. The order parameter then 
takes the form A(r) = A(r)e me . n = corresponds to a 
vortex-free state, and n= 1 to a vortex state. u(r) and f(r) can 
be expanded in terms of <j>j i. The coherence length may be 
estimated £(T) = Hvf/itA(T) [8], with v? the Fermi velocity 
and A the order parameter of the vortex-free state at R — > oo. 

In Fig. 1 (a) and (b), we plot the spatially averaged SC order 
parameters A. In Fig. 1 (a), for the vortex-free state (n = 0) at 
T = 0, the superconductivity remains robust as R decreases. 
The sudden drop in A at R <C £(0) is due to the quantum 
size effect, where the energy level spacing due to confinement 
becomes comparable to the SC gap [9]. The vortex state (n = 
1) vanishes at R < R c ps 1.5£(0), consistent with the GL 
theory [5] and the recent experiment [10]. Therefore, there is 
a sharp crossover between vortex and vortex-less behaviors at 
i?«1.5£(0). 






FIG. 1. (color online) s-wave superconductivity in a disk of radius 
R. (a) Spatially averaged SC gap A vs. R at T = 0. (b) A and the 
coherence length £ ({)) vs. T for R = 2.5£(0). ■ is for the vortex 
free state, and • for the vortex state in a magnetic field, (c) Phase 
diagram H vs. T for R = 2.5£(0). Here $o is the superconducting 
flux quantum. The parameters are gm/2-rrh 2 =0.256, fx= 180 meV, 
A = 30 meV and m = m c , the free electron mass, which give the SC 
gap A = 1.2 meV, and £(0) ~ 40 nm in the large R limit. 

At small R> R c , the vortex state exists at low T, but van- 
ishes at T = T' < T c , due to the increasing vortex core size 
(or the value of £(T)) with T, as shown in Fig. 1 (b). The 
corresponding phase diagram in Fig. 1 (c) can be calculated 
by comparing the Gibbs free energy Q = (H) — TS + Th 
of different states in equilibrium [11], with H the mean-field 

Hamiltonian, S the entropy and Th = J dr J — 87r 1 the mag- 
netic field exclusion energy [12]. There is a sharp drop for the 
upper critical field at T ps 0.6T C , indicating the suppression 



of the vortex state. This result may explain the recent exper- 
iment on ultra small aluminum thin films, showing a similar 
drop below T c in the H — T phase diagram [13]. 

We now study the much more intriguing p-wave case. We 
shall only consider states which are allowed by rotational in- 
variance and shall restrict our attention to the chiral p x + ip y 
phase, the 2D analog of 3 He-A [14]. For a given spin compo- 
nent, the BdG equations have been derived by Matsumoto et 
al. [15] and read 
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with II(r) = E± [A±D± + |(n±A±)] . Here A± are 
the pairings of p x ± ip y states, respectively, and satisfy the 
self-consistent equations 
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with fc F = y/2mn/H 2 and D± = e ±l0 {d r ± [16]. In 
the thermodynamic limit, there are two degenerate eigenstates 
with p x ± ip y pairings. In the finite system, the two pairing 
symmetries are mixed near the boundaries. We consider here 
the state where the p x + ip y pairing is the dominant compo- 
nent. In the disk geometry, the pairing parameters have the 
forms, A+(r) = A + (r)e' m0 and A_(r) = A_(r)e l ("+ 2 ) , 
where n is the vorticity. We focus on three cases: the vortex- 
free state (n = 0), negative-vortex state (n = —1), and 
positive-vortex state (n = 1). Let us consider the spin triplet 
case with s z = 0. Our results are shown in Fig. 2. For the 
vortex-free state, the pairing is dominated by the p x + ip y 
component A+, except near the edge where p x — ip y com- 
ponent A_ also becomes substantial. At the edge r=R, 
|A + | ps |A_|. The right panels in Fig. 2 show the density 
distributions of the supercurrents. Note that at the edge, the 
current density shows rapid oscillations with a wave vector 
2kp. In Fig. 2 (b), near the edge within the length scale £(0) 
[18], the supercurrent flows counterclockwise, while next to 
the edge current within the penetration depth, about 2£(0), 
a weak screening current flowing clockwise is driven by the 
Meissner effect. We also confirm that the bound edge states 
contribute to the former and the scattering states contribute to 
the latter [19]. 

Due to the broken time reversal symmetry, there are two 
types of vortices with n = ±1, depending on the direction 
of the applied magnetic field. Note that there is an induced 
supercurrent near the vortex core, whose direction is the same 
as that of the intrinsic chiral edge currents for n = 1 and is 
opposite to that for n = — 1 . 

Now we turn to the study of the quantum size effect due to 
finite confinement at T = 0. Contrary to the s-wave SC state, 
there is no vortex-free p- wave SC state at R < R c ps 1.5£(0), 
and the negative-vortex state in the p-wave state is robust 
against small size, as we can see from Fig. 3 (a). The super- 
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paring the Gibbs free energy Q as usual. It is clear that this 
phase diagram and its s-wave counterpart in Fig. 1 (c) form a 
sharp contrast. It shows both re-entrant and induced p-wave 
SC phases in various regions, both of which are absent in the 
s-wave counterpart. Specifically, between 0.3T C and 0.6T C , 
there is a re-entrant SC phenomenon as H increases from zero 
since it evolves between two different SC states by crossing a 
normal phase; however, the SC(n=-l) state becomes induced 
from 0.6T C to T c due to the absence of vortex-free state. Sim- 
ilarly, a phase diagram for R = 1.4£(0) < R c is presented in 
Fig. 3 (d), showing the complete disappearance of the vortex- 
free state and the re-entrant phenomenon. Thus in this case, 
field-cooled samples with H ~ $ /ttR 2 will exhibit super- 
conductivity whereas zero-field cooled samples do not. This 
prediction can be readily tested by using S^RuC^ microcrys- 
tals [20]. Although strong evidence suggests that the SC state 
of this material has odd parity [21], the null result on the ob- 
servation of the edge currents [22] seems to shed doubt on its 
chiral p-wave symmetry [14]. Therefore, it will be exciting to 
confirm our prediction of the induced superconductivity in a 
magnetic field, since it can provide a very strong evidence of 
the chiral p-wave SC pairing for Sr2Ru04. 



FIG. 2. (Color online) Spin triplet p x ± ip y state (s z = 0) in a disk 
of radius R = 12£(0). Left: Spatial variation of the order param- 
eters A+ and A_ (imaginary part). Right: Spatial dependence of 
the supercurrent density je(r). Upper panels (a), (b): vortex-free 
state; middle panels (c), (d): negative-vortex state; and lower panels 
(e), (f): positive-vortex state. The insets on the right panels illustrate 
the corresponding current flowing directions at the edge and near the 
vortex core. The parameters are gm/2-Kh 2 = 0.2, fj, = A = 16.32 
meV and m = m c , which give T c =1.4 K, and f(T = 0) « 70 nm, 
close to those of Sr2RuQ4 [17]. 



conductivity in the vortex-free state is destroyed by the super- 
current induced by the chiral motion of the Cooper pairs at the 
edge, which becomes dominant in a sufficiently small sample. 
In fact, the boundary region in the vortex-free p-wave SC state 
is similar to the vortex core in the s-wave state where the su- 
perconductivity is suppressed by the supercurrent around the 
vortex core. In the negative-vortex state, the supercurrents 
around the vortex and around boundary flow in the opposite 
directions (see Fig. 2 (d)), and partially cancel each other, 
which makes the superconductivity stable. This situation is 
similar to that of a vortex-antivortex pair [4]. 

To propose experiments to test our findings, below we study 
the T-dependence of SC states for the system of R = 2.4£(0). 
As can be seen in Fig. 3 (b), the T-dependence of the order 
parameters for different SC states show qualitative discrep- 
ancy. In the absence of any external magnetic field, the chiral 
p-wave superconductivity disappears above 0.6T C . This is be- 
cause that £(T) increases with temperature. If a magnetic field 
is applied to the system to induce a negative-vortex, the SC 
condensate may revive between 0.6T C and T c . To further il- 
lustrate the whole physical picture quantitatively, we can con- 
struct the corresponding phase diagram in Fig. 3 (c) by com- 
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FI G. 3. (Color o nline) Spatially averaged order parameter A = 
•y/|A+| 2 + |A_ | 2 as functions of disk radius R at T=0 (panel a), and 
as functions of T for a given R = 2.4£(0) (panel b), in a spin triplet 
p x ± ipy-wave superconductor with s z — 0. ■ is for the vortex-free 
state, and • is for the negative-vortex state in a magnetic field. The 
parameters are the same as in Fig. 2. Panels (c) and (d) are phase dia- 
grams for R = 2.4{(0) > R c and R = 1.4£(0) < R c , respectively. 
In both cases, the upper critical fields are close to $o/27r£(T) 2 . 

To understand the finite size effects in both the s-wave and 
p-wave pairings from the topological point of view, we can de- 
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fine a winding number associated with the SC order parameter 

as 

2tti Jc A 

where A refers to the off-diagonal part of the BdG equations 
and C is a circle around the origin counterclockwise. For the 
conventional s-wave superconductors, the winding number W 
is just the vorticity n. However, for the p-wave pairing with 
Px + ipy as the dominant component, we find that W = n + 1, 
i.e., the sum of the vorticity and the additional +1 for the in- 
trinsic chirality. Therefore, forthep-wave vortex-free state the 
winding number is 1, similar to the vortex s-wave state. These 
two states with W = 1 are found to vanish below a critical size 
comparable to the coherence length. On the contrary, both the 
s-wave vortex-free state and the p-wave negative-vortex state 
with W = survive even below the critical size [23]. These 
results are in good agreement with an intuitive analysis on the 
GL free-energy density. We only focus on the dominant gra- 
dient term |X>A| 2 [8], where the covariant derivative T> reads 

T> = -iV - 2A = -id r e r - (i-d e + 2A 9 )e e (9) 

r 

in a rotationally invariant system, and we set h=e=c=l. Here 
^dgA dominates and makes the SC state disfavored, as the 
size is reduced. To avoid this energetically unfavorable term, 
dgA has to vanish, i.e., W — 0, leading to the survival of the 
s-wave vortex-free state and the p-wave negative-vortex state. 

The above result may be extended to discuss the equal-spin 
pairing (ESP) phase in a spin triplet p-wave superconductor, 
which may support half-quantum vortices [24]. In this phase, 
there are two weakly interacting condensates with Cooper- 
pair spin configurations | ff ) and | H) coupled by the elec- 
tromagnetic field. The Cooper pair wave function is given by 
[2] 

*(r) cx (A t+ (r) e ±4fl |n> + A±+(r)|W» * ( Px + i Py ) 

+ e 2ie (A t _(r)e ±4e |n> + A;_(r)|W» * (p* ~ iPy), 

(10) 

where A a denotes the pairing in the state \<j<j) and * is the 
symmetrized product. Then the two condensates are described 
by two separated sets of self-consistent equations with differ- 
ent vorticities but the same vector potential. In this situation, 
electrons with up spin form a vortex state with = ±1, 
while those with down spin form a vortex-free state. Based on 
the relation between the winding number and the finite size 
confinement discussed previously, it is straightforward to pre- 
dict that when the size of the system is reduced below R c , the 
pairing in | H) breaks down, while the pairing in | yf ) with 
iVf = — 1 is still robust. In this sense, a spinless chiral p-wave 
superconductor is achieved. We have confirmed this scenario 
by performing numerical simulations. We stress that the re- 
sults above are based on the assumption that Cooper pairs with 
distinct spins are weakly coupled. 



In summary, we have used BdG equations to study the fi- 
nite size effect, for both s-wave and chiral p-wave supercon- 
ductors. For the p-wave pairing, the vortex-free SC state does 
not exist below a critical size, whereas the vortex state is ro- 
bust even for the system size as small as the coherence length, 
where the opposite winding of the vortex compensates the p- 
wave intrinsic winding at the boundary. These results pre- 
dict a magnetic-field-induced superconductivity in ultra small 
samples with p-wave pairing such as Sr2Ru04. For all these 
quantum geometrical constrains for both s-wave and p-wave 
pairings, the winding number W plays a determining role. Al- 
though our mean-field theory has neglected SC fluctuations, 
we expect that our results should be valid for BCS supercon- 
ductors such as Sr2Ru04 with a large coherence length and 
superfluid density. 
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